Eigenstate Randomization Hypothesis: 
Why Does the Long-Time Average Equal the Microcanonical Average? 
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We derive an upper bound on the difference between the long-time average and the microcanonical 
ensemble average of observables in isolated quantum systems. We propose, numerically verify, and 
analytically support a new hypothesis, eigenstate randomization hypothesis (ERH), which implies 
that in the energy eigenbasis the diagonal elements of observables fluctuate randomly. We show that 
ERH includes eigenstate thermalization hypothesis (ETH) and makes the aforementioned bound 
vanishingly small. Moreover, ERH is applicable to integrable systems for which ETH breaks down. 
We argue that the range of the validity of ERH determines that of the microcanonical description. 

PACS numbers: 05.30.-d, 03.65.-w 



in 



i 

C 

O 

o 



> 

(N 
cn 



X 



Statistical mechanics is among the most successful 
fields of physics and has an extremely diverse range of 
applications. However, why statistical mechanics works 
so well has only recently begun to be understood. To 
answer this question, we must address the problems of 
equilibration and thermalization. In other words, does a 
given time-dependent physical quantity relax to a certain 
constant value? If so, can the equilibrated value be cal- 
culated by invoking statistical ensembles? If a system in- 
teracts with energy-exchanging reservoirs, the canonical 
ensemble is known to be useful for describing the system. 
Several works 1-3] have exemplified that entanglement 
between the system and reservoirs plays a critical role in 
justifying the canonical ensemble. 

When the system is isolated, or if there is no reservoir, 
the microcanonical ensemble is useful for describing the 
system. Reimann [J] addressed the issue of equilibration 
and showed that time-dependent expectation values of 
observables take on their long-time averages (LTA) most 
of the time. However, to justify the method of the micro- 
canonical ensemble, we have yet to address the problem 
of thermalization, and understand why LTA is equal to 
the microcanonical ensemble average (MEA) . 

The problem of thermalization has been addressed 
mainly by two approaches. One is based on eigenstate 
thermalization hypothesis (ETH) 0-0] j which states that 
diagonal elements of observables in the eigenenergy ba- 
sis take on a constant value in the thermodynamic limit. 
If ETH holds, then LTA=MEA is satisfied. However, 
no proof of ETH has been given and numerical tests are 
hampered due to computational complexity. Moreover, 
ETH is known to break down in integrable systems 
The other approach is to assume a virtual probability 
distribution [8]. This approach assumes that the initial 
state of the system is generated according to the virtual 
probability distribution. In this case, LTA=MEA holds 
with a high probability. However, in reality, the initial 
state is not generated in such a manner, and therefore 



the physical ground of this approach is unclear. 

In this Letter, we show LTA=MEA by resolving the 
problem of thermalization. Firstly, we show an inequal- 
ity, which gives an upper bound on the difference be- 
tween LTA and MEA. Secondly, as a sufficient condition 
for LTA=MEA, we derive an eigenstate randomization 
hypothesis (ERH), which implies that the diagonal ele- 
ments of observables behave randomly in the eigenenergy 
basis. We show numerically that ERH holds for few-body 
observables whether or not ETH holds. Moreover, we 
show that the microcanonical description is justified not 
just for expectation values but also for higher moments 
and the extent to which ERH holds delimits the range 
of validity beyond which the microcanonical description 
breaks down. Lastly, we show analytically that ERH 
holds for almost all observables. 

We consider an isolated many-body quantum system, 
which is described by Hamiltonian H, and a few-body 
observable A. Here by few, we mean a number much less 
than the number of particles in the system. Let the eigen- 
values and eigenvectors of H be E\ < E 2 < ■ ■ ■ < £jv , 
and \Ei) , \E 2 ) , . . . , \En ), respectively. We assume that 
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if and only if Ei t — Ej 1 and Ei 2 = Ej 2 . 



This condi- 



tion is necessary for the many-body quantum system to 
equilibrate according to the theorem found by Reimann 
and its generalization @, 0|. If H has some degeneracies, 
we choose {\Ei)}^ 1 to diagonalize A in each degenerate 
subspace. The initial state is expanded in terms of 
energy eigenstates as — J^i^i 1^)- We assume that 
the coefficients {Ci} i= ?j satisfy the following condition. 
There exists a macroscopically small energy width AE 
such that d ^ only if Ei 6 [E - AE, E + AE], where 
E = (t/j\H\tp). We assume the number of eigenenergies 
in the interval [E — AE, E + AE] is much greater than 
one, yet AE is much smaller than the resolution limit of 
an experimental apparatus. 
Now, LTA of the observable A is defined as LTA^ (A) = 
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limT^oo T" 1 / Q (ip(t)\A\ip(t)} dt. Based on the assump- 
tions described above, we can show that 

LTA^A) =]T|Q| 2 ^, (1) 

i 

where Aj = (Ei\A\Ei). Given an energy window 8 (> 
AE) [loj], MEA of the observable A is defined as 

MEA S (A)=N(E,5)- 1 ]T A h (2) 

Bi-_E|<5 

where AT (E, 8) is the number of energy levels that satisfy 
\Ei — E\ < 8 and will be denoted simply as N. Note that 
N can be as large as N ~ 10 10 for macroscopic systems. 
For the sake of convenience, we relabel i so that {Ei}^^ 1 
appear in [E — 8, E + (5] in ascending order. As shown 
later, the inequality 

|LTA V ,(A) - MEAj(A)| < a A D ( ^ + ^ , (3) 

which gives an upper bound on the difference between 
LTA and MEA, holds for any positive integer to (< N). 
Here, a\ is the variance of the sequence {A,}^ 1 , a\ = 

N- 1 Eto 1 ^ ~ ^) 2 > where A = A" 1 4, °i m Y 

is the variance of coarse-grained A, which is defined below 
Eqs. (@| and ([5]), and D is a measure of the smoothness 
of the wave function such that | \C i+1 \ 2 -\C l \ 2 \< D/N 2 . 

Now, we define our coarse-graining procedure and 
cr(m). Given a positive integer to, we divide N by to 
and define the quotient M and the remainder n as fol- 
lows: M = \N/m~\, or M is the smallest integer not 
less than N/m, and n = N — (M — l)m. We group the 
elements of the sequence {A;}^^ 1 to by to and define 
another sequence {-Bfcj^Lo 1 : 

^ m— 1 

B fe = -VA fcm+Q (fc = 0,l,...,M-2), (4) 

TO ' — ' 

n— 1 

%-l = — ^ A( M -l)m+a- (5) 

We call {Bfelfcl^ 1 the m-th coarse-grained sequence of 
{Aij^Q 1 . The variance of the new sequence is ct(to) 2 = 
M-^m^ -S) 2 , where B = M' 1 ^^ B k- 
Then, the first coarse-grained sequence is {A;}^ 1 it- 
self, so cr(l) = cta. Note that if the sequence {Ai}^ 1 
behaves randomly, fluctuations of its m-th coarse-grained 
sequence are greatly suppressed and its variance a(m) 2 
is expected to decrease as to increases. 

Let us consider the meaning of the inequality (|3|), 
which consists of the summation of two terms: a(m)l a a 
and m/N. As we make further coarse graining and in- 
crease m, the first term <j(ra)j(JA decreases depending on 
the randomness of the matrix elements of the observable 
A. On the other hand, the second term increases and 



its magnitude is determined by the smoothness of the 
weights of the wave function. So, in making the upper 
bound of ([3]) minimal, we have a trade-off between the 
randomness of the observable and the smoothness of the 
wave function. 

Now, let us consider sufficient conditions for 
LTA=MEA. If ETH, which asserts that the first fac- 
tor a a vanishes in the thermodynamic limit, holds, then 
LTA=MEA holds. However, there exists another suf- 
ficient condition that a(m) decreases rapidly as to in- 
creases. So, we come to ERH, which asserts ct(to) = 
tr^TO -7 , where 7 is a positive parameter. As explained 
before, this implies that the sequence {Ai}^ 1 behaves 
randomly. Note that, if ETH holds, ERH automatically 
holds, since <t(to) = for every to. In this sense, ERH 
includes ETH as a special case. If ERH holds, by min- 
imizing the upper bound of ([3]) with respect to to, we 
get 

|LTA^(A) - MEA 5 (A)| < <x A D'N-^, (6) 

where D' = 7 1 /(7+i)(i + ^^D. If the weights {\C l \ 2 } l 
are smooth and D = O(N ), the right-hand side is 
vanishingly small. Note that the requirement on the 
smoothness of weights is much weaker and actually D — 

(AT-7/(7+l)). 

To show numerical evidence for the behavior ct(to) = 
ct^to -7 (7 > 0), wc adopt a hard-core Bose-Hubbard 
model: 

(.id) (*d) 

where stands for the nearest neighbors. The oper- 
ators b\ and bi are the creation and annihilation opera- 
tors of site i, respectively, and satisfy [bi, bj] = [&],&]•] = 
[b h b)} = for i ^ j, {bi,bi} = {bl,bt} = 0, and 

{bi, b]} — 1. The operator rii = b\bi is the number opera- 
tor of site i. In the following we measure energy in units 
of J and take U = 0.1. We take 21 sites and arrange 
them on a 3 by 7 rectanglar lattice. We denote the lat- 
tice constant as d, and each lattice point as Xi — (l x d, l y d) 
(l x = 0, 1, . . . , 6 and l y = 0, 1, 2). Since the Hamiltonian 
conserves the total number of particles N p , the dynam- 
ics proceeds in the subspace with a definite value of N p . 
Here, we choose A p = 5 due to the limitation of com- 
putational complexity. The dimension of the subspace is 
( 2 5 x ) - 2 x 10 4 . These settings are to be compared with 
those by Rigol et al. 0, who studied ETH. 

The observable A we choose here is the quasi momen- 
tum distribution 

n(k x ,k y ) = -±-J2 eiiXi ~ Xi)Pb t b j> (8) 

* 1,3 

where L x and L y are the sides of the rectangle (L x = Id 
and L y = 3d) andp = 2ir{k x / L x , k y /L y ) (k x = 0, 1, ... ,6 
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FIG. 1: (Color online) Energy eigenvalue Ei (in the 
units of hopping parameter J) versus the correspond- 
ing eigenstate expectation value of the observable n(0, 0) 
{{Ei\n{0,0)\Ei)). Globally (Ei\n(0,0)\Ei) varies smoothly 
(upper figure), whereas locally it oscillates randomly around 
its local median (lower figure). This local behavior reflects 
the fact that even if two adjacent energy eigenvalues are very 
close, the corresponding many-body wave functions are quite 
different. 



and k y = 0, 1,2) is the wave vector. Here, we have set 
h = 1. The average energy and energy width of the initial 
state are chosen as E = and AE = 0.1. The energy 
width is chosen so that the central value of the sequence 
{^4i} i=u does not vary in the region [E — AE, E + AE]. 

The relation between energy eigenvalue Ei and the cor- 
responding eigenstate expectation value Ai is illustrated 
in the lower figure of Fig. [TJ If we focus on a small energy 
width, Ai behaves randomly around a certain value as ex- 
pected. The decay of the variance a(m) = UAmT 1 caused 
by this random behavior is illustrated for four choices of 
n{k x , k y ) in Fig. [2j where 7 is found to be between 0.41 
and 0.50. 

Changing the configuration of the system from the 
rectangle to the one-dimensional loop, we have performed 
the same analysis. By this change, the system becomes 
integrable 11]. In this case lattice points are denoted 
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FIG. 2: (Color online) The variance of the m-th coarse- 
grained sequence a(m) 2 plotted for four momentum distribu- 
tions n(k x ,k y ) for a non-integrable system: (k x ,k y ) = (0,0) 
(cross), (0,1) (asterisk), (1,0) (filled square), (1,1) (open 
square). The moment a(m) decays as er(ra) = <r(l)m -7 . The 
exponent 7 of the four observables are 7 = 0.50, 0.41, 0.45 
and 0.49, respectively. The decay in the moment occurs be- 
cause the original sequence (lower figure of Fig. [TJ behaves 
very randomly. 
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FIG. 3: (Color online) The variance of the m-th coarse- 
grained sequence a(m) 2 plotted for four momentum distri- 
butions n(m) in the integrable system: k = (cross), 1 (as- 
terisk), 2 (filled square), 3 (open square). The exponent 7 is 
found to be 7 = 0.26 (k = 0), 0.25 (k = 1), 0.26 (Jfc = 2) and 
0.31 (k = 3). 



as Xi = Ud (Ji = 0,1,..., 20) and the quasi momentum 
distribution is 



(9) 



where p = 2irk/L and L = 20d. The result is also illus- 
trated for four choices of n(k) in Fig. [3l where we have 
chosen E = and AE = 0.2. Although ETH is known to 
break down in this system Q, <r(m) decays as in the non- 
integrable case. Therefore, ERH holds and LTA=MEA 
is satisfied. In this case we find 7 to be between 0.25 
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and 0.31. Since {Ai}^^ 1 oscillates very drastically, ETH 
breaks down, but ERH holds because of the rapid oscilla- 
tions. Thus, even if ETH breaks down, LTA=MEA still 
holds. By these analyses, we find numerically that the 
moment decays as a(m) — aAin^' 1 (-y > 0) whether the 
system is integrable or not. 

We have numerically showed that ERH holds for few- 
body observables. However, ERH does not hold for all 
observables. The projection operator onto the z-th en- 
ergy eigenvector Pj = \E/) (Ei\ violates ERH but this is 
a full TVp-body observable. Let us assume ERH holds 
for less than JT max -body observables. Let A be a K- 
body observable, then A n (n — 1,2,...) is an nK- 
body observable and LTA^(A n ) - MEk AE {A n ) (n 

1.2 [Kmax/ K\). Our results demonstrate that the 

microcanonical description is valid not only for expecta- 
tion values but also for such higher moments. 

In addition to the numerical evidences, we can show 
ERH holds for almost all observables. Precisely speak- 
ing, if we consider the uniform measure on the set of 
N x N Hermitian matrices whose trace norm is gN 2 , 
A = {A TrL4U] = gN 2 ,A^ = A}, ERH holds for al- 
most all matrices, where g is a constant of the order of 
unity and the factor N 2 implies that every element of 
A does not vanish or diverge as N increases. The set of 
matrices A can be considered as an (N 2 — l)-dimensional 
sphere whose radius is y/ gN 2 . Note that (j4|) and §5§ 
are linear transformations such that they are written in 
vector representation as B = C(m)A, where C(m) is 
the M x N matrix. The variance cr(m) 2 is calculated 
to be cr(rn) 2 = M^ 1 A T C(m) T (I - M^ 1 \m)C (m)A = 
A T W(m)A, where 1m is the M x M matrix whose 
elements are all equal to unity. Then, (a(m) 2 ) = 
gTr[W(m)}, and (ct(to) 4 } = 5 2 (l+2/Y 2 )- 1 {Tr[W/(m)] 2 + 
2Tr[VL(m) 2 ]}, where (...) means that the average is to be 
taken over the (TV 2 — l)-dimensional sphere. For the sake 
of simplicity, here we consider the case in which N is di- 
visible by m, or n = m and N = Mm. The eigenvalues of 
W(m) are and 1/N, and their degrees of degeneracy are 
N — N/rn + 1 and N/m — 1, respectively. When m <C N, 
(a(m) 2 ) ~ gm" 1 and y / {a(m) 4 ) - (a(m) 2 ) 2 / '(cr(m) 2 ) ~ 
y / 2m/N < 1. These show that ERH holds for almost all 
observable A with 7 = 1/2. 

Finally, we prove (J3J). We denote g, = \Ci\ 2 , 
Then, the m-th coarse-grained sequences of {qi} 1 ^ 1 , 
{rk}k=o'i is defined as r = C(m)q. First, we use 
the triangular inequality and Schwarz's inequality to 
show \LTA^{A) - MEA AE (A)\ < J2 k a AA km+a (r k - 
km-\-a (qkm+ a - rk)\ < a(m)[(N + 
m)(Q(m) - 1/N)} 1 / 2 + a A [N(Q - Q(m))}^ 2 , where 
Q(m) = J2ka r 1- According to Cauchy-Schwarz's in- 
equality, 1/N < Q(m) < Q(l). Moreover, by summing 
both sides of |<& —Qj\ < \i — j\D/N 2 over i and j, we have 
2(AQ(1) - 1) < D 2 . From \q km+a - r k \ < Dm/2N 2 , 
Q - Q{m) < D 2 m 2 /N 3 . This proves ©. 



In conclusion, we have proposed a new scenario (ERH 
with dH)) under which LTA=MEA holds in isolated 
many-body quantum systems. Together with the work 
by Reimann [ij and its generalization this scenario 
explains why the method of the microcanonical ensemble 
gives the correct description of quantum many-body sys- 
tems. Our scenario is valid whether or not ETH holds. 
We verify ERH numerically in two typical cases of the 
hard-core Bose-Hubbard model. It is expected that ERH 
holds for other configurations and observables of this 
model. In fact, it is suggested that ERH holds in chaotic 
systems [l2|,[l3j]. Moreover, we have analytically shown 
that ERH holds for almost all observables. Once we find 
that ERH holds for up to if max -body observables, we 
also know that it sets the range of validity of the micro- 
canonical description and this knowledge can be used as 
an operational criterion of macroscopic observables. 
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